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Abstract. A new method, AMORE - based on a genetic algorithm optimizer, is presented for the automated study 
of colour-magnitude diagrams. The method combines several stellar population synthesis tools developed in the 
last decade by or in collaboration with the Padova group. Our method is able to recover, within the uncertainties, 
the parameters - distance, extinction, age, metallicity, index of a power-law initial mass function and the index 
of an exponential star formation rate - from a reference synthetic stellar population. No a priori information is 
inserted to recover the parameters, which is done simultaneously and not one at a time. Examples are given to 
demonstrate and to better understand biases in the results, if one of the input parameters is deliberately set fixed 
to a non-optimum value. 
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The appendix and the tables (electronically available at 
the CDS) listed therein are published in the electronic 
version of this paper. 



al. |1999[ Holtzman et al. |1997| & |1999|) the full, unbinned 



1. Introduction 

The present data flow of many ongoing surveys 
such as 2MASS (Beichma n et a l |1998[ Skrutskie |1999| ), 
DeNIS (Epchtein 
da Costa [1997 



OGLE-II (U dalski 
SDSS (Fan [1999 



1999), EIS (Rcnzinifc 



et al. [L9971 

da Costa 
et al. 

and references cited therein), and 



1997, da Costa et al. 
1997, Paczyhski et al. 



199S) 



1999) 



1998 



even upcoming surv eys as GAIA (Gilmore et al 
Perryman et al. [2001 ) — is so large that one requires cither 
a semi-automated or a fully automated method to analyse 
the colour-magnitude diagrams (CMDs) in the resulting 
databases. In this paper we discuss the development and 
the tests of an automated analysis method, which fully 
employs the colour and magnitude information available 
about the stars populating the CMD. Our method is based 
on an i mplem entation of the CMD diagnostics suggested 
by Ng ( 1998 ). The method uses, in contrast to other tech- 



niques (see Bertel li ct al. |1992| , Gallart et al. |l996| fc |199S 



Geha et al. 1998 Harris & Zaritsky 2001, Hernandez et 
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distribution of magnitudes and colours of the stars popu- 
lating the CMD. 

The purpose of this paper is to verify that astrophys- 
ical parameters for a synthetic single stellar population 
can be reliably retrieved with the so-called 'AutoMatic 
Observation REnderer' AMORE. In Sect. [2] an outline of 
AMORE is given together with its individual building 
blocks. In Sect. |^ we outline the method we use and in 
Sect. ^ we describe the tests performed with synthetic stel- 
lar populations. The results are given in Sect. || and we 
discuss in Sect. || the practical limits on the convergence, 
which is imposed by some degeneracy of the parameter 
space. We end with prospects on forthcoming tests, recom- 
mendations for improvements, and an outlook on future 
developments. 



2. AMORE 

2.1. Project outline 

AMORE tries to find the best matching synthetic CMD 
to an observed CMD. Such a synthetic CMD contains 
for stellar aggregates the contribution of one or more 
stellar populations at the same distance. In the case of a 
CMD along a particular line of sight in our Galaxy the 
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synthetic CMD can moreover contain the contribution 
of various populations with stars distributed at different 
distances. 

In this paper we focus on the implementation and the 
performance of AMORE for the fitting of CMDs. For 
sake of argument only one, single stellar population has 
been considered. The implementation of automatic fitting 
CMDs with multiple stellar populations with stars at the 
same or different distances will be subject of forthcoming 
papers. 



2.2. Building blocks 

AMORE, see Fig. [l] and Sect. for details, combines 
various analysis tools developed and improved at Padova 
during the last decade. Conceptually, it is made up out of 
the following building blocks: 



— a synthetic Hertzsprung-Russcll Diagram gen erator 
(hereafter referred to as HRD-ZVAR see Sect. 2.3), 

— a mockup version of the HRD Galactic Software 
Telescope (HRD-GST, see Sect. j2~4[ ), 

— a statistical diagnostic tool to compare t he ob served 
and synthetic CMDs with each other (Ng 199§|) , and 
the PIKAIA (extended v ersi on) and POWELL optimizers 

which search for 



(see respectively Sects. f2J3| and |2 



the best fit between an observed and a synthetic CMD 
in a multi-parameter space. 

In the following subsections a description is given of the 
various building blocks. 

2.3. HRD-ZVAR 

In the late-eighties Bertelli developed a code to generate 
synthetic Hertzsprung-Russcll diagrams (HRDs) from the 
isochrones computed by the Padova group (cf. Chiosi et al. 
1989] ). Initially the synthetic colour-magnitude diagram 



(CMD) technique was applied mainly in the studies of 
LMC open clusters^ (see for example Bertelli et al. 1985, 
1990| or Chiosi et al. |1989 ), through which the amount of 



convcctivc overshoot was calibrated for the computation of 
a new generation of stellar evolutionary tracks. Successive 
improvements were gradually applied when new sets of 
evolutionary tracks (see Bertelli et al. 1994 for details) 
were compu ted with improved radiative opacities (Iglesias 
et al. |1992| ). 

The backbone of HRD-ZVAR, the extended version of 
the HRD generator, is formed by the evolutionary tracks 



computed by Bertelli et al. (199C; Z = 0.001), Bressan et al 



1 The analysis technique with stellar ratios was employed. 
The reason for this is that ratios are less sensitive to uncer- 
tainties in certain regions in the CMD, which might not be 
reproduced properly due to various reasons such as the in- 
put physics used for the calculations of the stellar evolutionary 
tracks or the transformations from the theoretical to the ob- 
servational plane. 



( 1199$ Z = 0.020), and Fagotto et al. (1994abc; Z = 0.0004, 
0.004, 0.008, 0.050, 0.10). The mctallicities of the tracks 
follow the enrichment law AY/AZ = 2.5 (see references 
cited in Chiosi 1996 and Pagel & Portinari 1998). 



HRD-ZVAR indicates that the mctallicity Z is not lim- 
ited to the fixed values for which the evolutionary tracks 
have been computed, but is variable through interpolation 
between the metal-poorest and metal-richest tracks avail- 
able inside the database of evolutionary tracks. In this way 
one is able to generate synthetic stellar populations with a 
smooth metallicity coverage. A prerequisite however is to 
use a complete and homogeneous library of evolutionary 
tracks and some improvements are expected if one adopts 
a grid of tracks with a smoother metallicity coverage. 

HRD-ZVAR has been distributed (privately) to various 
research groups and is also referred to as ZVAR|. The 
version d istribu ted, mod ified, and used by f or example 
Aparicio (|l999| ), Gallart ( |l998|) , and Ng et al. (|l996| , |l997|) 
is from now on referred to as VI. 0. Version VI. 6 is used 
for the simulations and results presented in this paper. 
This version contains a number of modifications and im- 
provements which speed up the code and fix some (rarely 
encountered) bugs which interfered with the automatic 
minimization process. Although various analysis methods 
are available, we limit ourselves here to the description 
of the parameters related to the HRD-ZVAR as adopted for 
AMORE. 

After selection of a set of tracks with fixed metallicity 
and the choice of the parameters t]rgb and ijagb, (the 
mass loss along the Red Giant Branch (RGB) and the 
Thermally Pulsing Asymptotic Giant Branch (TP-AGB) 
phases respectively), the major input to be specified for 
HRD-ZVAR are: 

— the metallicity range, Z ranges from Z m j n to 7i rnax ; 

— the age range, the age ranges from t m j„ to t max ', 

— the slope a for the power-law IMF (Initial Mass 
Function); and 

— the index (3 for the exponential SFR (Star Formation 
Rate) . 

The luminosity and effective temperature for each syn- 
thetic star of arbitrary metallicity is transformed to an 
absolute magnitude in a photometric passband with the 
metho d out lined by Bertelli et al. ( 1994 ) and Bressan 
et al. ( |1994| ). Default setup for HRD-ZVAR is the UBVRI 
JHKLMN^j broadband photometric system. The setup can 
be altered to mimic any system, given the description of 

2 Note that an older and modified version of HRD-ZVAR is 



actually the program used to generate the Bertelli et al. (1994) 
isochrones. HRD-ZV AR is free of the interpolation difficulties as 
reported by Olsen ( |l999| ) 

3 The SFR is only fixed for the time being to the adopted 
exponential shape 

4 Note that the IR photometric system is based on an 'aver- 
age' photometric system as described by Bessell & Brett (1988). 



Proper transformations ought to be applied to the actual pho- 
tometric system prior to any astrophysical interpretation of the 
results. 
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Fig. 1. AMORE flowchart: schematic diagram of 
the individual building blocks. PIKAIA outlines 
the direction of the evolution paths to be in- 
vestigated. Input for the stellar population syn- 
thesis engine HRD-ZVAR is the Padova library of 
evolutionary tracks. The luminosity and effec- 
tive temperature for each synthetic star of arbi- 
trary metallicity is then transformed to an abso- 
lute magnitude in a photometric passband with 
the method outlined 
and Brcssan et al 



by Bertelli et al. fll994| ) 
The synthetic HRD is 



( |1994| ). 

then 'observed' and 'detected' through a Monte 
Carlo 'observing run' with the HRD-GST. Note 
that an alternative route is possible with single 
stellar populations (SSPs). The synthetic CMD 
is compared with the observed diagram and a fit- 
ness parameter is subsequently communicated to 
PIKAIA, which suggests a new set of parameters 
for each trial. The iteration lasts for a user de- 
fined, fixed number of trials. POWELL's method 
of minimization is subsequently applied to get 
closer to the local or global minimum. After com- 
putation of the uncertainties for each parameter 
the evolutionary run is either aborted or a new 
PIKAIA cycle is started after shrinking the limits 
of the parameter space (see Sect. 2.9) 



the spectral response of the filter and the detector of the 
photometric system. 



2.4. HRD-GST 



HRD-ZVAR was integrated in a galactic model by Ng (1994, 
1997ab). The HRD-GST (Galactic Software Telescope) has 
been applied in the studies of the galac tic st r uctur e to- 
wards the G alactic Cen tre (Ng et al. 1995 , 1996 and 



ll. 


1995 


1997 





Bertelli et al. 1995 , 1996] ) and other regions in our Galaxy 



(Ng et al. |1997|) . In this paper we do not require the 
full complexity of the structural properties from the GST 
model. We only use a limited number of options to 'ob- 
serve' a synthetic HRD at the suggested distance and to 
simulate the photometric errors, extinction and crowding. 

A table of the photometric errors, covering a specific 
magnitude interval per passband, is used and the program 
interpolates linearly to obtain the intermediate values. We 
assume for the simulation that the photometric errors are 
Gaussian distributed. A different description of the pho- 
tometric errors will be used when published artificial star 
tests (Stetson and Harris 198S , Gallart et al. 1999] ) on an 
observational data set are indicative for a significant de- 
viation from a Gaussian behaviour. 

The visual extinction is simulated with the average value 
provided and appropriately scaled to a value in different 
passbands. In the UBVRI passbands we adopted the scal- 



ing according to van dc Hulst (1949; curve no. 15) and for 



the JHKLMN passbands we follow the scaling laws pro- 



vided by Rieke & Lcbofsky ( 1985 ). We do allow for some 
random scatter around the average extinction. However, 
we do not consider (yet) the effects due to patchi ness o f 
the extinction along the line of sight. Ng& Bertelli ( 1996 ) 
demonstrated that this is in first approximation, visually 
almost indistinguishable from a random scatter around an 
average extinction. 

In many studies the observations are crowding limited, 
due to the increasing number of stars towards fainter mag- 
nitudes. Crowding gives rise to star blends which affects 
the magnitude and the colour of the stars. The group of 
stars will be detected as a single star with a magnitude 
equal to the sum of the stellar flux of the stars involved in 
the stellar blend. The remaining stars are 'hidden' from 
detection. 

The blends are well described as unresolved, apparent bi- 
naries. The simulation of apparent binaries is made with 
an itcratively improved blending probability, which is de- 
fined as the probability that a star within a given ensemble 
of stars might blend with another star from the same pop- 
ulation. Each synthetic star within a stellar population is 
tested against the blending probability. 
The percentage of artificial binaries is with this scheme 
about twice the blending factor. The blending factor in 
different passbands is not necessarily the same and the 
occurrence of star blends is furthermore not necessarily 
correlated, due to possible differences in the exposure time 
or seeing conditions. 

The fainter companion stars of artificial binaries will give 
rise to incompleteness of the synthetic stellar sample. This 
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allows us to map the synthetic stellar completeness func- 
tion, which can be compared directly with the complete- 
ness function obtained from artificial star tests from the 
observed stellar sample. 

In the CMD to be analysed we assume implicitly that, 
between the observed and synthetic photometric system, 
the uncertainty in the magnitude zeropoint is smaller than 
the uncertainty in the zeropoint of the colour, see Carraro 
et al. (1999), and references cited therein. We allow for this 
reason the possibility that a small zeropoint shift might 
be present between the colours of these systems. 

2.5. SSPs 

Although straightforward, the generation of a large num- 
ber of synthetic HRDs with each their own specific age- 
metallicity range, star formation history and initial mass 
function can be a time consuming task, because of the 
repetition of many calculations to generate one diagram. 
Figure [l] indicates that an alternative route with single 
stellar populations (SSPs) is available for the automated 
analysis. However, this method requires the computation 
of a large, regular grid of SSPs for different age-mctallicity 
ranges and IMFs. The star formation history and age- 
mctallicity range are the result of the linear combination of 
the SSPs. The time spent on computations of new CMDs 
thus can be greatly reduced through the use of proba- 
bility density diagrams. However, this is not our prime 
objective. Our present goal is to develop an automated 
fitting method, which compares on a star by star basis, 
and to demonstrate its potential. Optimization for speed 
is not yet our primary concern. Moreover, the generation 
of one synthetic CMD with 5000 stars takes about 0.75s 
on a PC equipped with a 200 MHz Pentium processor. 
This is a good indication that our present version of the 
software tool is performing at an acceptable speed. We re- 
fer to Dolphin ( p97[ |200l| , p002| ) and Olsen ( p999[ ) for a 
description of a method using SSPs. 



2.6. PIKAIA 

The PIKAIA optimization packa ge (V 1.0; public domain) 
was developed by Charbonneau (1995) and a full descrip- 
tion of this package is given by Charbonneau & Knapp 
( |1996| ). PIKAIA has been used successfully in a wide 



range of astrophysical applications (e.g. Bobingcr 2000 
Charbonneau et al 



Kaastra et al 
et al 



1998. Gibson & Charbonneau 



1996, Kcnnclly et al. 



1996, Mcintosh et al. 1998, Metcalfe 1999 



1998 



1996, Lamontagnc 
Mcwe 



et al. [1996| , Noyes et al. |1997| , Saha |199S| , Wahdc [1998|) 



PIKAIA is based on a genetic algorithm (Holland 1975, 
Goldberg |j~989| , Davis |1991| , De Jong |1993D , and is in prin- 
ciple not a function optimizer, but it does this extremely 
well. It searches for, locks on to, and pins down an opti- 
mal solution in a way, which is conceptually comparable to 
biological evolution through natural selection. Genetic al- 
gorithms arc capable to explore and find in a robust way 



an optimum, but not necessarily the best, setting for a 
particular problem. In our case this comes down to mini- 
mizing the difference between a synthetic and an observed 
CMD by evolving the astrophysical parameters that define 
the shape of the CMD (see Table | and Sects |lo| and |). 
We follow the generally accepted biological terminology 
for the description of a genetic algorithm. 

A genetic algorithm makes use of a reduced version 
of the evolutionary process. The gene pool, i.e. the set of 
parameters to be optimized, and its associated phenotypic 
population evolves in response to 

— the reproductive success of the population, following 
its reproduction plan; 

— genetic recombination (crossover, see Fig. ||) at breed- 
ing; 

— random mutations (see Fig. ^J) during breeding, which 
affect a subset of the events. 



2.6.1. Flow control 

PIKAIA has 12 flow control parameters which arc discussed 



by Charbonneau (1995) to whom we refer for a detailed 
description. We limit ourselves to a short summary and 
the purposes of these control parameters: np defines the 
number of individuals in a population; ngen specifies the 
number of generations that the population is evolving; nd 
is number of digits encoding accuracy^ used for the pa- 
rameters; pcross is the probability that a crossover oc- 
curs between the chromosomes of the parents; imut, pmut, 
pmutmn, and pmutmx specify the mutation mode, the mu- 
tation rate (the initial mutation rate if the rate spans the 
range from pmutmn to pmutmx for imut = 2), and the mini- 
mum and maximum mutation rate; f dif the fitness differ- 
ential controls the selection of the individuals for breeding 
through their fitness; irep defines the reproduction plan 
to be followed; ielite defines if the fittest individual can 
or cannot be selected for replacement; and ivrb specifies 
verbose mode for extra on screen information during the 
evolutionary run. 

Table [l] holds a list of the PIKAIA flow parameters which 
were kept constant during all the simulations described in 
this paper. In the following subsections we describe the 
extensions added to the 12 flow parameters of PIKAIA. 

2.6.2. Crossover 

The crossover operator is very effective in a global explo- 
ration of the full parameter space and is in a way compara- 
ble to a variational calculus method. A one-point crossover 
scheme, see Fig. is sometimes inadequate to combine 
and pass on c ertain features encoded on the chromosomes 
(Michalewicz 1996 ) to its offspring. In some cases only 
a correlated modification of a number of genes, say 2, 



5 The 1 or 2 number of digits encoding accuracy maps the 
parameter via an integer value, either 0-9 or 0-99, to a range 
of floating point values by a controlling 'normalization' func- 
tion, see Charbonneau (1995) for additional details 
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Fig. 2. Schematic view of the main genetic operators acting on two sets of parameter strings in order to procreate 
two new sets (one chromosome set for each individual) of test strings: natural selection of two individuals G and <j 
for the breeding of offspring, modification of the genetic content by means of a homologous crossover - at gene Gi/ gi, 
and finally the mutation of gene g% to *<?2- Each gene on a chromosome represent a parameter to be optimized. The 
schematics for a 2-point crossover scheme are shown for comparison 



Table 1. PIKAIA control parameters. These parameters were kept fixed at the values listed throughout the evolutionary 
runs 



initial 


PIKAIA 


parameter 


comment 


value 


default 


identifier 




too 


too 


np 


the number of individuals per population 


20 


500 


ngen 


number of generations 


2 


5 


nd 


number of digits encoding accuracy 


2 


2 


imut 


mutation mode, imode=2 then pmut= [pmutmn. pmutmx] 


0.005 


0.005 


pmut 


initial mutation rate 


0.005 


0.005 


pmutmn 


minimum mutation rate 


0.35 


0.25 


pmutmx 


maximum mutation rate 


0.95 


1.0 


fdif 


fitness differential 


3 


1 


irep 


reproduction plan 





1 


ielite 


elitism 








ivrb 


verbose mode 



will result in a fitter offspring. This can, for example, be 
mimicked through the application of a two- or multi-point 
crossover scheme. PIKAIA has been extended with the con- 
trol parameter rcross to handle a multi-point crossover 
operation. For example: rcross = 1 represents the default 
one-point crossover, while rcross = 2. 3 represents a '2.3- 
point' crossover: i.e. a two-point and three-point crossover 



for respectively 70% and 30% of the cases (2.3 = 0.7x2 + 
0.3x3). 

2.6.3. Brood recombination 

An important drawback of genetic algorithms is that the 
crossover operator is for about 75% of the time lethal to its 
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offspring, i.e. it produces children which arc not as fit as 
their parents (Banzhaf ct al. 1998| ) . To avoid missing, po- 
tentially, fitter offspring and to reduce the destructive ef- 
fect of the crossover operator we have incorporated brood 
recombination in PIKAIA through a new control param- 
eter rbrood. The default PIKAIA reproduction scheme is 
obtained with rbrood = 1, i.e. two parents breed once and 
produce two new individuals. With rbrood = 3. 5 the par- 
ents are allowed to breed on average 3.5 times to produce a 
larger offspring (in this case 7 on average). For rbrood > 1 
one ends up with more than two offspring. PIKAIA on the 
other hand expects from each pair of parents only two 
children. This constraint was obeyed in order not to sig- 
nificantly alter the global behaviour of PIKAIA. Therefore, 
in order to avoid an exponential growth of the population 
for rbrood^ 1, only the fittest two individuals of the local 
offspring survive^ and enter the global population for a fit- 
ness evaluation from which a selection is made for further 
breeding. 

The extra breeding increases the computational effort con- 
siderably, due to a larger number of function evaluations. 
The advantage is a more rapid increase of fitter individu- 
als through the selection of effective crossovers from good 
recombinations. On the other hand, a rapid increase of fit- 
ter individuals might lead to a premature convergence to 
a local minimum, due to a smaller variance in the genetic 
pool. 

2.6.4. Creep mutations 

We further intro duced the 'creep' mutation (Charbonneau 
&Knapp 1996 ) in order to overcome the so-called 
'Hamming Wall' problem, i.e. the inability to cross in a 
decimal encoding scheme certain boundaries with a one- 
point mutation operator. The creep parameter pcreep de- 
fines the probability that a gene in the pool undergoes a 
'standard' mutation (change digit randomly in the range 
-9) or the 'creep' mutation (add or subtract one from the 
current value of the digit). Wc adopted as default equal 
weight for the occurrence of a 'creep' or 'standard' muta- 
tion. 



2.6.5. Correlated mutations 

In general mutations occur in the optimization process 
to avoid premature convergence. A low mutation rate is 
sufficient for this purpose. However, a high mutation rate 
can be used as an additional way to explore the parame- 
ter space like a virus. Although one would prefer to use a 
(multi-point) crossover operator, we do allow that muta- 
tions can be used instead, simply because the two opera- 
tors also co-exist in nature. 



This is irrespective if the local, fittest individuals are 
weaker than the weakest individual in the global population 
or if the locally remaining offspring are fitter than the fittest 
individual in the global population 



We modified PIKAIA's uniform mutation mode. In the 
majority of the cases we require a (anti-)corrclatcd change 
between two or more parameters (see Sect. |6.2[ ). In a stan- 
dard mutation scheme convergence might be slow if one 
has to wait for the simultaneous occurrence of a favourable 
(anti-) correlated mutation of two specific parameters in 
order to improve the fitness. We introduced an extra pa- 
rameter pcorr which defines the probability that a corre- 
lated mutation occurs. If this is not the case the standard 
mutation scheme is chosen. Otherwise we allowed that in 
50% of the cases the mutations of two genes (igenl and 
igen2 are extra input parameters added to the modi- 
fied version of PIKAIA) are more relevant than the mu- 
tations occurring in other parameters. For the remaining 
50% of the cases the two genes are determined stochasti- 
cally. Additional details about the adopted values of the 



control parameters are given in Sect. 2.1C and Tabic A.l 



2.7. Fitness 

PIKAIA searches for the optimum solution by maximizing 
the fitness function /. The fittest solution has a fitness 
/ = !, while the worst has / = 0. We use the Ng (F99Sh fit- 



ness function, a combination of a chi-squared and Poisson 
like functions. These functions minimize the differences 
between observed and synthetic diagrams via a star-by- 
star matching^] scheme. 
The Ng fitness function is defined as: 



/ = 



1 



1 + F 
where F is 



F = Fl + Fl 



(1) 



(2) 



F x is the chi-squared function of the best fitted points 
within a3~5<7 error ellipse and Fp is the Poisson function 
of the residual points outside this ellipse. Both F x and 
Fp are dimensionless, but they hold information about 
the average uncertainty in units of cr, say a x and ap. For 
example: the average uncertainty per point for a fit with 
F x is F x x a x . 

F x and Fp are respectively defined as: 



F x = Vx 2 = Vx 2 (0,S)/N match , 



and 



F P = 



No,not + Ns^not 



(3) 



(4) 



The intuitive motive behind this is to make a division be- 
tween the synthetic points matching the observed CMD 
for which the errors are expected to be normally dis- 
tributed and the points which do not match and are 
allegedly assigned to the Poisson merit function. The 
method actually uses F x as a loosely fixed 'anchor', puts 



This option is feasible due to the increment of the present 
day computational speed 
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Table 2. AMORE control parameters, used as the 'educated guess' for the astrophysical parameters in the first PIKAIA 
cycle. These values result in a fitness / = 0.004. The value for the parameters arc constrained between the lower and 
upper limit 



# 


parameter 


value 


lower limit 


upper limit 


description 


1 


log d 


3.7 


3.6 


4.0 


log of distance 


2 


A v 


0.30 


-0.001 


0.5 


extinction 


3 


log t iow 


9.5 


9.0 


10.3 


log lower age limit 


4 


[Z]ioiu 


-0.90 


-1.69897 


0.69897 


log lower metallicity limit 


5 


a 


1.35 


1.001 


3.5 


IMF slope 


6 


log thigh 


9.6 


9.0 


10.3 


log upper age limit 


7 


[Z]hgh 


0.30 


-1.69897 


0.69897 


log upper metallicity limit 


8 


P 


1.0 


-2.00 


5.00 


SFR slope 



the outlier points in Fp and then reduces the number 
of unm atched points by minimizing Fp. We refer to Ng 
( |l998|) for additions details and a discussion of these func- 
tions. Suffice to say that for an acceptable solution both 
F x and Fp are about 1 or smaller^ which on its turn can be 
relaxed to the condition F^ 2 [] and thus yields ^ ^ / < 1. 

The formal lcifc uncertainty of each parameter k, see 
Table ^, is obtained through variation of this parame- 
ter and by minimizing the function | \/Fk — V F m i n — 1| . 
Conceptually, this is similar to moving the merit function 
in the Fp , F x — plane away from its optimum setting, to 
the nearest position on a contour +luk higher. The asso- 
ciated fitness function /o-.fe^l is: 



fcr.i 



II 



(5) 



where k is the particular parameter for which the uncer- 
tainty is estimated and F min is the global value obtained 
for the fittest population. 

2.8. POWELL 

We implemented a hybrid optimizer in which we use 
PIKAIA to explore the parameter space and then use 



POWELL's minimization algorithm (Powell 1964; Press et al 



1986) to pin down the nearest local or global minimum 



through a direction set method which produces N mutu- 
ally conjugate (non-interfering) directions. For details and 
an excellent description of this algorithm we refer to Press 



8 For our testcase, as described in Sect. ^, F x ranges from 
0.7-1.0 for all observed data points with a matching synthetic 
point within a 3a uncertainty ellipse. This corresponds to a 
goodness of fit parameter ranging from 0.49-1.0. It further 
indicates that it is justified to assume that the measurement 
errors are normally distributed. 

9 This condition is comparable b ut no t equivalent to the re- 
sults obtained by Gallart et al. (1999). They demonstrated 



from a comparison with colour-magnitude bins that a good 
agreement between the input and recovered SFR(t) required a 
reduced chi-squared of xt — 2-0. 

10 For both age and metallicity ranges the associated uncer- 
tainties denote for the lower and upper values the — la and 
+ la boundary 



et al. and references cited therein. 

A hybrid minimization strategy is used, because PIKAIA 
is by definition not a function optimizer, but it tends to 
get close near a (local) optimum. POWELL is used to get 
even closer to the (local) optimum. If we had landed in 
a local optimum then we needed PIKAIA to jump out of 
it. The origin of our need for a hybrid search strategy is 
comparable to the minim ization problems encountered by 
Harris & Zaritsky J200l[). 



2.9. Contracting parameter space 

A full exploration of parameter space at one digit accuracy 
would take considerable time (in the order of weeks), even 
though certain forbidden combinations of parameters can 
be excluded d priori. A full exploration at two or more 
digits accuracy is nearly impossible due to present day 
computational limits. 

We implemented a dynamic, scalable parameter range 
in our search for an optimum set of parameters. The pa- 
rameter range shrinks after each optimization cycle with 
POWELL. This leads to an improved accuracy in the results 
with a fixed number of digits encoding accuracy. 
The automated re-scaling of the parameter range improves 
the resolution of the exploration of the search grid. In ad- 
dition, due to the re-scaling one may circumvent partial 
degeneracy of the parameters. 

The global function F(F x ,Fp) gives the global dis- 
tance to the minimum in terms of a 2 (u Xl ap). If we have n 
parameters then each parameter k is in a simple approach 
on average about \/F/n o~k away from its optimum value, 
because we assume that 



F 



fe=i 



Vk.tr 



(0) 



where yk.sim is the simulated value of parameter k and 
Uk,true is the true value of parameter k. The average off- 
set per parameter k is therefore J F/n. The new limits 
can then be set to i-^/F/nov To balance the cancel- 
lation of errors, due to a negative correlation between 
some parameters, we adopt a rather conservative ap- 
proach by constraining the limits of each parameter k to 
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±3 x 1.3^/F/nak- In addition we add the condition that 
3xl.3V^7™>0.005. 

2.10. Adjustable parameters: free &. fixed 

AMORE uses in essence two different sets of parameters 
as input. One set contains program flow control parame- 
ters and information about the observational data to be 
simulated: the photometric errors & crowding factors per 
passband, a shift of the zeropoint of the colours due to a 
difference between the observed and synthetic photomet- 
ric system, and the spread around the average extinction 
value to account partly for the differential reddening in a 
field. 

The other set of input parameters is used by PIKAIA 
and can be divided into two parts (see Table [I] and 
Table ^). One part contains the PIKAIA control parame- 
ters. The other part contains the lower and upper limits of 
the astrophysical parameters to be optimized as well as an 
initial guess for the value of those parameters (first column 
of Table resulting HRD in Fig. ||b). These parameters 
are a combination of the synthetic population's intrinsic 
properties, i.e. age, metallicity, slope of the power-law IMF 
and the index of the exponential SFR. In our case we al- 
low for age and metallicity not to be restricted to one fixed 
value, but to cover a specific range (see Table ||). 
In addition, there are two parameters which mimic the 
synthetic population's behaviour as placed in a mockup 
version of our Galaxy. These parameters arc the distance 
from the Sun and the average extinction. 

In total we thus have 8 free parameters which AMORE 
has to optimize simultaneously. 



3. Method 

3.1. Genetic evolution 

In terms of genetic programming the objective of AMORE 
is to determine the genome (i.e. the set of astrophysical 
parameters described in Sect 2.10 , see also Table ^) of a 
specified individual (i.e. the observed stellar population). 
Note that it is not possible to directly observe the genome. 
The genome is determined from the phenotype of each in- 
dividual (i.e. the synthetic CMDs, sec for example Fig. ^). 
The genetic information is located in the genes of one 
chromosome]^] 

A guess of the genotype of the observed CMD is 
obtained through comparison with a synthetic CMD, 
which is generated via supervised evolution and breeding 
(PIKAIA together with POWELL). The stars in the synthetic 
CMD population with a particular genotype are raised to 
maturity (HRD-ZVAR and HRD-GST). A group of individu- 



11 The genetic information is currently located on one chro- 
mosome. Individuals with two chromosomes might be consid- 
ered as a future extension. A two-chromosome approach has the 
advantage that certain genetic information can remain present 
in a recessive form. 



als^is allowed to procreate (the chance of an individual 
procreating depends on its fitness and the selection pres- 
sure, see Charbonneau 1995| and Charbonneau & Knapp 
1996) and the genetic information of the parents is passed 
on to their offspring (see Fig. ||). 

A fitness evaluation (a comparison between the ob- 
served and synthetic CMD) provides a ranking of the re- 
sulting group of individuals. If the individual has 'good 
genes' it survives, remains in the group and gets a chance 
of procreation. 

The evolutionary process of breeding and fitness eval- 
uation is repeated for a fixed number of generations. The 
gene pool of the resulting best individual at the end of 
the evolutionary run with AMORE hopefully represents a 
near-optimum representation of the unknown genome. 

3.2. Running AMORE 

Initially PIKAIA is in control (see Fig. [l]) of the evolu- 
tion for a fixed number of generations. Afterwards POWELL 
tries to improve the genome of the fittest individual com- 
municated through PIKAIA. We then determine the un- 
certainty for each gene on the chromosome. Subsequently, 
we tighten the limits on the range of variation allowed for 
each gene and re-scale the parameters on the genetic print 
of the fittest individual accordingly. During the shrinkage 
of the parameter range we do not re-scale the genetic in- 
formation of the remaining individuals, but preserve their 
former values as semi-random input for the continued op- 
timization process. The latter addition to the hybrid 
scheme is most likely a significant driver in speeding up 
the search for a fitter individual. 

After each optimization with POWELL a new cycle with 
PIKAIA is started with the current best parameter set as 
'educated next guess' for AMORE's progressive evolution. 
The tot al n umber of PIKAIA cycles is user defined. 
In Sect. |2.9| w e argued that the parameters are on average 
about y/F/no~k away from its optimum value. The con- 
vergence however is not governed by the average 'distance' 
that each parameters is away from its optimum setting. It 
is mainly determined from the ability to tune the param- 
eter which has the largest offset from its optimum value. 
In the AMORE training sessions it was noted that with 
F ~3.0 about three of the eight parameters are about 1 
(~ J F/3(Jk) away from their optimum value. AMORE 
has a built-in option to do a random variation from - 3 o~k 
of two parameters (randomly selected) from the running 
best solution when the fitness is less than 0.30. Above this 
threshold we choose a new value for two parameters ac- 
cording F J pcrfc, where p can be any number between 
1.0 and 4.0. 



4. Tests 



12 To avoid confusion the term group is used instead of the bi- 
ological term population, because each individual in the group 
is actually a synthetic stellar population 
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1.5 



0.5 <?: 



1 2 



Fig. 3. Evolutionary statu s for all the trial models 
(filled dots; see Table A.l for details) for each pop- 
ulation after 400 generations, see also Fig. || The 
filled star -fc points to the location obtained with 
the original parameter settings as given in Table |^. 
Acceptable solutions are found in the region for 



which both F x and Fp are less than 1, see Sect. [2.7| . 
The condition has been relaxed to F < 2 (solid white 
line). The asterisk indicates the position obtained 
for a different realisation of the test population by 
changing the random seed. The shaded regions indi- 
cate solutions for which the difference between the 
CMDs from the 'observed' and synthetic population 
is on average less than respectively 3 *7, y3 a, and la. 
Note that rounding errors, see Sects. |t|&U, give 
rise to a degeneration of the parameter space around 
/>0.25 (i.e. F<3 ) 
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Fig. 4. The filled circles in panels a-f display the values obtained for the parameters using the models in Table A.l 
The solid line refers to the value set for the original population 



4.1. Test objectives 

We performed several tests on AMORE in order to 

— verify and validate AMORE's performance in retriev- 
ing the astrophysical parameters of a synthetic single 
stellar population; 



— determine adequate values for the parameters pcross , 
rcross, rbrood, pcreep, and pcorr in the ex- 
tended version of PIKAIA; 

— study the effects of rounding and degeneracy; 

— study the effects of fixing parameters on the conver- 
gence; 
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— study the effects of a high extinction on the conver- 
gence. 

4.2. Setup 

The hybrid interaction between PIKAIA and POWELL, com- 
bined with a progressive shrinking/re-scaling of the pa- 
rameter space, requires that a trade-off has to be made in 
the choice of the size of the population and the number of 
generations we allow this population to evolve in order to 
obtain results in a reasonable amount of CPU processing 
time. 

We explored several different settings for the PIKAIA con- 
trol parameters, because the tuning of those parameters is 
very problem dependent (Charbonneau & Knapp, 1996| ). 
The values we decided to use are listed in Table |l. Four 
notes can be made here. 

Firstly the steady-state-delete-worst reproduction plan 
(irep=3) we adopted, in which we replace the least-fit in- 
dividual from the population when the fitness of the new 
individual is superior to that of the least-fit population 
member. Choosing this reproduction plan implies that the 
elitism control parameter (ielite) is non-operative, be- 
cause elitism is active by default. We evaluate d two other 
reproduction plans (Charbonneau & Knapp 1996 ); full 
generational replacement and steady-state-delete-random. 
The steady-state-delete- worst reproduction plan produced 
on average the best results. 

Secondly the mutation rate of 0.35 corresponds, in case of 
a default 2 digit accuracy, with the on average occurrence 
of 2.8 mutations per astrophysical parameter. 
Thirdly, the fitness differential parameter fdif , a mea- 
sure for the selection pressure, would normally be chosen 
as high as possible (f dif =1 in this case). However, it may 
possible to circumvent local minima by lowering that value 
a bit (Charbonneau & Knapp, 199£). Setting fdif =0.95 



turned out to be a good trade-off choice. 
Fourthly, we want to explore as large a fraction of the pa- 
rameter space as possible at the first entry in AMORE. 
This is done by using only an one digit accuracy (nd=l). 
Due to the active re-scaling of the parameter space bound- 
aries we do not require a very high precision in our explo- 
ration. A one percent accuracy (nd=2) of the parameter 
space is sufficient in the subsequent PIKAIA cycles. 

In biological terms, the PIKAIA control parameters de- 
fine the ecosystem in which our population evolves. 
All computations presented in this paper were performed 
with an executable generated with the g77 compiler. This 
executable was then installed on various PCs running Red 
Hat Linux0 6.X and 7.X. The PCs were equipped with 
Intel Pentium III or Athlon processors with clock speeds 
ranging from 600-1200 MHz. 

All tests, unless stated oth erwis e, use the synthetic 



population as described by Ng ( 1998 ): 
— a metallicity range, spanning Z = 0.005 -0.030; 



13 Red Hat® is a registered trademark of Red Hat Software, 
Inc. and Linux® is a registered trademark of Linus Torvalds. 



— an age range from 8-9 Gyr; 

— an initial mass function with a Salpeter slope; 

— an exponential decreasing star formation rate with a 
characteristic time scale of 1 Gyr. 

The test population contains N = 5000 stars and is placed 
at 8 kpc distance. The 'observational' limits are set to 
V U m=22 m and l lim = 21 m . 

4.3. Description of the tests 

4.3.1. Test 1: Determining values for pcross, 
rcross, rbrood, pcreep, and pcorr 

In the fir st t est, we evaluate the 162 models listed 
in Table [A. l| in order to study the effect of the 
PIKAIA parameters pcross , rcross, rbrood, pcreep, 
and pcorr on the convergence and computational effort. 
The test has as a secondary objective to provide an un- 
derstanding of the degeneracy of the parameter space. 
All astrophysical parameters to be retrieved are set free, 
floating between reasonable minimum and maximum val- 
ues (see Table H for details). AMORE runs for 20 iterations 
of 20 generations (ngen=20) to recover the a priori known 
parameters of the synthetic population. The number of 
iterations and generations determine the total length of 
an evolutionary run: 20x20 = 400 generations. Note that 
the range of each parameter is set within reasonable lim- 
its and not taken excessively large, because it might lead 
to the case that no acceptable parameter setting is found 
with the standard iteration loop. 

4.3.2. Test 2: Rounding 

The second test deals with the effects of rounding. We vary 
the number of significant digits in the input parameters 
to reveal AMORE's sensitivity to rounding. In this case 
we do not make an evolutionary run, because the trial 
set of parameters (the educated guess at the start-up of a 
PIKAIA cycle) is the correct one. For clarity we label the 
fitness in this test by f trial instead of Ja, the fitness after 
a complete evolutionary run. 

4.3.3. Test 3: Fixing parameters at the correct value 

In the third test we take six models in which one of the pa- 
rameters is set fixed at its correct value in order to study 
the effects on the convergence. The models chosen were 
two of high, two of intermediate and two of low fitness 
as determined from the first test. The convergence in this 
test basically can go two ways: either the convergence is 
faster, because less parameters have to be optimized. Or, 
due to the fact that AMORE has less maneuverability in 
this situation, the convergence is slower. We adjusted the 
limits for age and metallicity as given in Table|| such that 
AMORE would not try to find solutions in forbidden re- 
gions of parameter space which might severely slow down 
convergence due to constant rejection by AMORE of the 
chosen parameter values. 
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Table 3. Average fitness (/a) for dif feren t values of the 
parameters as obtained from Table A.l . The first col- 
umn displays the parameter name which is varied, the 
parameter value is given in the second column, the third 
Sz fourth column show respectively /a, uncorrected for 
pcorr=0.0, together with its standard deviation, and 
the fifth and sixth column display the averaged values af- 
ter removal of the results of the models with pcorr=0.0. 
See Sect. 15.1.2 



for additional details 



parameter 


value 


fA 




fA 


Cn-l 


pcross 


0.50 


0.276 


0.070 


0.287 


0.066 




0.85 


0.299 


0.063 


0.297 


0.056 


rcross 


1.00 


0.293 


0.065 


0.290 


0.063 




2.00 


0.285 


0.068 


0.297 


0.065 




3.00 


0.285 


0.070 


0.290 


0.057 


rbrood 


1.00 


0.286 


0.077 


0.230 


0.059 




2.00 


0.292 


0.058 


0.301 


0.055 




4.00 


0.285 


0.066 


0.276 


0.067 


pcreep 


0.0 


0.287 


0.071 


0.292 


0.068 




0.3 


0.290 


0.065 


0.297 


0.058 




0.7 


0.285 


0.066 


0.288 


0.060 


pcorr 


0.0 


0.278 


0.077 


NA 


NA 




0.3 


0.285 


0.065 


NA 


NA 




0.7 


0.300 


0.057 


NA 


NA 



Table 4. The effect of rounding in the parameters on the fitness evaluation; f trial is the fitness obtained with AMORE 



for the given values of the parameters. The average value and a are based on the 162 models in Table A.l. The a is 
obtained from an unweighted average for the fitness of the models involved. The a is the error as estimated by Ng 



199S 



parameter 


log d(pc) 


Ay 


logti ow 


log t hgh 




[Z}hgh 


a 


P 


f trial 








original 


3.906335 


o m oo 


9.90309 


9.95424 


-0.60206 


0.17609 


2.35 


1.0 


0.44597 


round- vl 
round-v2 
round-v3 


3.906 
3.906 
3.9063 


o m oo 

0™00 

o ra oo 


9.903 
9.903 
9.9031 


9.954 
9.954 
9.9542 


-0.60 
-0.602 
-0.602 


0.18 
0.176 
0.176 


2.35 
2.35 
2.35 


1.0 
1.0 
1.0 


0.28595 
0.30812 
0.42439 


average value 

a 

a 


3.8958 
0.0033 
0.012 


m 027 
m 014 
m 06 


9.866 
0.049 
0.043 


9.984 
0.047 
0.023 


-0.554 
0.050 
0.18 


0.244 
0.13 
0.08 


2.358 
0.034 
0.03 


1.574 
1.40 
1.4 





For example, fixing the Z,i ow parameter at its correct value 
of -0.60206 means that we have to adjust the lower limit 
for Z n igh to -0.60206 as well. 

In the case of fixing the log ti ow parameter this also im- 
plies that the initial guess has to be adjusted. We set this 
initial guess to 10.1. 



4.3.4. Test 4: Fixing parameters at the wrong value 

In the fourth test we take six models in which one of 
the parameters is set fixed at ler offset (determined from 
the first test) from its original value, in order to study 
its effect on the 'second best' setting of the remaining 
parameters. Normally one would expect a fitness />■§■ 
In this case, however, F< Fp + F 2 = l 2 + (1 + l) 2 = 5 
and the associated fitness constraint drops to / > g • 
However, this assessment ignores the fact that, when a 
parameter is offset from its optimum value, the num- 
ber of matched points will decrease and Fp increases. 
Using Eq. (^) one has for a good fit F = 2. On average 
the offset per parameter k from the optimum value is 

o~k = \ o~k , at best the offset is <7fc , and in the 

worst case this is s/2ak- So with one parameter k 
put at 1 <7fc offset we distinguish the three possibilities 

1 at best F=(l + 2) =3 -> /=| 

2 on average F=(l + ±) 2 + | =4 -► f=\ 

3 at worst F=(l + \/2) 2 = 5.8 -> / = 



Note that the worst case limit is in agreement with the 



results presented in Table A. 3 



The effect of the 1 a offset of one of the parameters will 
partly be canceled by forcing other parameters away from 
the optimum value. For example, the effect of an increased 
extinction can be masked partially by generating a bluer 
stellar population with a lower metallicity and a younger 
age. The effect will be such that the fitness will not be 
around / ~ , but somewhere in the range g^g < / < | . 
We fixed the parameters both at one sigma above and one 
sigma below the original value, because the evolutionary 
effects do not have to be symmetric. The only exception 
is the extinction, which we only fix at one sigma above 
the original value of Ay = 0.0 

Again we adjusted the limits for the upper and lower 
limit for age and metallicity. 

5. Results 

AMORE has been tested for a wide range of setups. The 



results in Table A.l indicate that AMORE give both ac- 



5 

l 

(S.N 



ceptable and less acceptable solutions. They are displayed 
in Figs. | & f§ 

To better understand what goes on during the ge- 
netic evolution we display the results from model [A.l[ -40. 
Figure || displays an example of the evolution of the merit 
function F for a number of generations. It shows how 
the initially dispersed individuals gradually find their way, 
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start to cluster together around generation 10, and pene- 
trate the region with acceptable solutions after about 50 
generations. After 100 generations the improvements be- 
come marginal for this model. 

Figure ^| displays for the same model |A.l[ -40 the phenotyp- 
ical changes of the CMD for several fitnesses during the ge- 
netic evolution. The various panels show that the synthetic 
CMD resembles better and better the 'observed' CMD 
when the fitness improves. Note that at fitness / = 0.05 
one already gets for the eye appealing solutions. 
Figure [?] shows the improvements of the astrophysical pa- 
rameters as a function of increasing fitness for the models 
[A.l| -40 and A. 1-51. The panels for distance and extinction 
show that the distance is systematically underestimated, 
while the extinction is overestimated. But in general one 
notice s tha t the astrophysical parameters obtained from 
model |A.l[ -40 get quite close to the parameters of the CMD 
to be matched. 



5.1. Test 1: Parameter values and degeneracy 
5.1.1. Degeneracy of the parameter space 



Table AT is displayed in Fig. |3j. The clustering in the 
figure provides an indication that a degeneracy of the pa- 
rameter space is present near / > 0.25 (i.e. F < 3, see 
cq. ([!])). Without a major computational effort it will be 
difficult to obtain a significant improvement of the param- 
eters once / > 0.25. However, F < 3 indicates a region in 
the (F x , Fp) -plane for which the systematic offset of the 
individual parameters from its true value are on average 
less than y/F/n tJk = 0.6erfc, see Sect. 2.9 for details. In 



practice it turns out that a strong correlation between 
three of the eight parameters has the culprit; at least two 
of them have to change simultaneously in the proper di- 
rection in order to improve the fitness (see also Sect. |6.3| ). 
They have an average offset of ~ \J F/3 au = 1 Cfc, while 
for the remaining parameters this is « 1 ak- 
in addition, Fig. |] displays the retrieved parameters 
for all models as a function of fitness. Note, that AMORE 
systematically underestimates the distance of the test pop- 
ulation. On the other hand, the effect of this underestima- 
tion is in its turn partially canceled by overestimating the 
extinction, the upper age limit and the slope of the power- 
law IMF slightly (see also Figure [j]) . Another clue we get 
from Fig. |^ is that the slope of the SFR (3 is very poorly 
constrained. 



5.1.2. Determining values for pcross , rcross, 
rbrood, pcreep, and pcorr 



Looking at Table A.l, a result that is immediately obvious 
is that pcorr = 0.0 has a strong stabilizing effect on the 
simulation. Its effects overwhelm the effects of the pcreep 
parameter and lock j 'a at a certain value, which may 
be quite good (e.g. models 10, 11 and 12) or quite poor 
(e.g. models 28, 29 and 30). Setting pcorr > is thus 
preferred to avoid getting locked in a low value for Ja- 



Table 5. Fitness statistics when fixing one parameter at 
its correct value. Averaged fitness values j a and their as- 
sociated standard deviation er rl _i are obtained from sim- 
ulations with the setup parameters from models 9, 14, 22, 
34, 40 and 52. See Table [0 for additional details 



parameter 


Sa 


Cn-l 


log d 


0.299 


0.081 


Ay 


0.274 


0.081 


logtlow 


0.332 


0.021 


log thigh 


0.304 


0.056 




0.269 


0.036 


Tjhigh 


0.361 


0.036 


a 


0.305 


0.066 


P 


0.312 


0.076 



In order to compensate for this strong stabilizing effect, 
we also evaluate in Table || the average fitness of the 
models when we exclude all models which have pcorr = 
0.0. 

As expected, the rbrood parameter has a strong influence 
on the amount of computational time needed. Although 
the models with high values of rbrood are somewhat 
better than models with low values, this effect is only 
marginal. Considering that a high value of rbrood lessens 
the genetic variation in the gene pool while increasing the 
computational time needed for a run with several factors, 
it is desirable to have a low value of rbrood. 



The different parameters are not independent, as 
can be seen from Table |A.l| a nd Table |[ Simply taking 
the best options in Table |3j yields model 134 for the 
case in which pcorr=0 . has not been corrected for, a 
reasonable, but not an exceptionally good model. 

5.2. Test 2: Rounding 

Table |] shows the effect of the accuracy of the retrieved 
values of the parameters on the evaluation of the fitness. 
If one applies a rounding to one or two significant digits it 
is not possible to reproduce the expected fitness, i.e. the 
expected fitness drops from f trial = 0.44 to f trial =0.28. A 
better agreement can be obtained by reporting the values 
of the parameters with the addition of one or more (ap- 
parently insignificant) digit (s). A renewed search (Table IJ 
round-v3) circumvents the local optimum at f tr i a i = 0.28 
and results in a near global fitness of ftriai = 0.42, which 
is close enough to the value expected. 
The true a line in the table shows that both the Z 



high 



and 



the /3 parameter are the weak links in the overall param- 
eter estimation (see also Fig ||) . 

5.3. Test 3: Fixing parameters at the correct value 



The results of fixing parameters at the correct value are 
listed in Table |A.2| and an example of the diagnostics is 
listed in Table la Details of the individual setups for these 
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Fig. 5. Conception diagram of the evolution of the genetic population of model [A.l[ -40 during the optimization 
process displayed in Fig. |?]. Frame a shows the initial population and the frames b-1 show the population after several 
generations up to generation=400. The outer shaded region indicates solutions for which the difference between the 
CMDs from the 'observed' and synthetic population is on average between l-3<7. The inner shaded regions marks 
the region with solutions for which the difference between the 'observed' and synthetic CMDs are less than 1 a. Such 
solutions are close to perfect matches between the 'observed' and synthetic data and are considered to belong to a 
group of solutions for which one may say "too good to be true". The solid lines indicate the 10 a, 20 a, 30a and 40 a 
contours 



tests are given below. In general, the results of the tests 
for which one of the parameters was set to the correct 
value were slightly better than the results for the models 
for which all parameters are set free, see Tables || & || for 
additional details. This behaviour is due to the fact that 
by forcing one parameter to a fixed value the evolutionary 
path changes. The models were selected from the results 



with low and intermediate fitness given in Table A.l 



5.3.1. Fixed distance and extinction 

Distance and extinction are negatively correlated. When 
both parameters are set free, a certain degree of degen- 
eracy is to be expected. Fixing one of the parameters at 
its correct value can break this degeneracy. The result de- 



pends strongly on the evolutionary path of the other pa- 
rameters. 

The lower value for the average fitness /a = 0.299 ± 0.081, 
when fixing the distance at the correct v alue , is caused 
by the presence of one outlier (see Table |A.2| ), which is 
caused by the age-metallicity degeneracy. Excluding this 
value results in an average fitness of /a = 0.328 ± 0.039. 
In general: the extinction can be reliably retrieved when 
fixing the distance. 

When considering a fixed extinction, the results show a 
strong variation in both age and metallicity. It should also 
be noted that the average 10 log (distance (pc)) retrieved 
is only 3.8953 ± 0.0066. This is more than one sigma away 
from the optimum value for the distance (see Table ||). 
This is an indication that retrieval of the distance by fixing 
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Fig. 6. Genetic evolution of the colour-magnitude diagram (CMD) from the first test population. Panel a displays the 
riginal population to be matched. The physical parameters for this population are described Table ^ and Sect. 4.2 



o 



The CMD of the initial trial population is shown in panel b. Panels c — e display the resulting CMDs obtained with 
setup |A.l[ -40 for different fitnesses (see Sect 2.7). The fitnesses / = 0.05, / = 0.19 and / = 0.27 are respectively reached 
after 20, 60 and 80 generations. The dots in the panels b — e are used for each matching point, while the red open stars 
in panels c — j are the points in the simulation which have no counterpart in the original CMD. Panel f displays the 
fitness / = 0.41 as obtained after 341 generations. Note that the CMDs of panels a and f as well as c — e are visually 
almost indistinguishable. Panels g— j displays the residuals between the simulated and the original CMD (a; green 
solid squares ■) are those points in the original which have no counterpart in the simulation 



the extinction is hampered by the agc-metallicity degen- 
eracy Therefore, the distance cannot be reliably retrieved 
when fixing the extinction to its correct value. 



5.3.2. Fixed age and metallicity 

Fixing one of the age limits results in values for both the 
age and metallicity which are close to the input values. 
This is due to the (partial) breaking of the age-metallicity 
degeneracy. The distance-extinction degeneracy remains. 
The results also suggest that the age-metallicity degener- 
acy has a stronger impact on the fitness than the distance- 
extinction degeneracy. 



Fixing the upper metallicity limit to its correct value 
shows that the values for age and metallicity come closer 
to their original, input values. This is quite in contrast 
with the results obtained from fixing the lower metallicity 
to its correct value. Table A. 2 shows that both the high 
metallicity limit and the slope of the exponential SFR are 
not well constrained. This behaviour can be accounted 
to the implicit shape of the linear age-metallicity relation 
adopted in the HRD-GST. The number of high metallicity 
stars is smaller than the number of low metallicity stars 
due to the adopted, exponentially decreasing ([3 = 1), star 
formation rate. The consequence is that the high metal- 
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Fig. 7. Panels a-h display the convergence curves for the parameters of models [A.l[ -40 (dotted line; fitness / = 0.41) 
and |A.l[ -51 (long dashed line; fitness / = 0.25). The solid line in the frames a — f refers to the value adopted for the 
original population. The long, dot dashed line in frames i&j shows the threshold values to be crossed for acceptable 
solutions, i.e. F < 2 and -| < / < 1. The short dashed area in frame j marks the region where degeneracy of the 



parameter space becomes noticeable (sec Sects. 6.3 for details) 



licity limit can be better determined when both the low 
& high limits are determined in union. 



5.3.3. Fixed slope for the power-law IMF 

Fixing the slope a of the power-law IMF to its correct 
value ensures stability in the magnitude direction of the 
CMD. This implies that the degeneracy in the colour de- 
pendent parameters, like age, metallicity, and partly the 
star formation rate, becomes more apparent. Although 
the overall results are quite good, only model |A.2[ -14 is 
significantly affected by this degeneracy. If we leave out 



model A. 2-14 from the statistics the average fitness be- 



comes f A = 0.331 ± 0.020. 



5.3.4. Fixed slope for the exponential SFR 

The exponential star formation rate parameter [3 is tied 



to the age & metallicity range, see also Sect. 5.3.2. Fixing 



the parameter (3 better constrains in particular the upper 
metallicity limit. However, it does not avoid that the ge- 
netic evolution enters into a local age-metallicity gap, see 



Table A. 2. Excluding model A. 2-34 improves the average 
fitness in Table | to /a = 0.342 ± 0.029. 



5.4. Test 4: Fixing parameters at the wrong value 

The results of fixing the parameters at a ± ler offset from 



its original value are listed in Table A. 3. An example of 
the diagnostics of these tests are given in Table |8|. Details 
of individual setups are given below. 
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Table 6. Description of the diagnostic statistics for model [A.l[ -40 when fixing one parameter at its correct value, see 



Sects. 5.3,2.7 and 6.2.1 for additional details 



model 


No, not 


Ns,not 


-^V match 




F P 


F 


fA 


ideal 










5000 


0.000 


0.000 


0.000 


1.000 


free 




60 


60 


4940 


0.847 


0.849 


1.438 


0.410 


fixed: 


log d (pc) 


155 


155 


4845 


0.903 


2.192 


5.621 


0.151 


fixed: 


A v 


103 


103 


4897 


0.840 


1.457 


2.827 


0.261 


fixed: 


logti oro (yr) 


76 


76 


4924 


0.817 


1.075 


1.823 


0.354 


fixed: 


logih 9 h(yr) 


104 


104 


4896 


0.902 


1.471 


2.978 


0.251 


fixed: 




94 


94 


4906 


0.836 


1.329 


2.466 


0.289 


fixed: 


\Ahgh 


76 


76 


4924 


0.815 


1.075 


1.820 


0.355 


fixed: 


a 


74 


74 


4926 


0.834 


1.047 


1.799 


0.357 


fixed: 


P 


80 


80 


4920 


0.727 


1.131 


1.808 


0.356 



Table 7. Fitness statistics when fixing one parameter 
± la from its original value. Averaged fitness values 
and their associated standard deviation <t„_i are obtained 
from simulations with the setup parameters from models 
9, 14, 22, 34, 40 and 52. See Table [A^ for additional de- 
tails 



parameter 


offset 


fA 


(7,i-l 


log d 


-la 


0.291 


0.048 




+ 1(7 


0.280 


0.066 


Ay 


+lcr 


0.284 


0.020 


\0g t low 


-1(7 


0.258 


0.044 




+ 1(7 


0.279 


0.012 


log thigh 


-1(7 


0.229 


0.022 




+ 1(7 


0.207 


0.014 




-1(7 


0.254 


0.066 




+ 1(7 


0.315 


0.027 


%high 


-1(7 


0.270 


0.054 




+ 1(7 


0.283 


0.018 


(X 


-1(7 


0.266 


0.025 




+ 1(7 


0.273 


0.027 




-1(7 


0.232 


0.071 




+ 1(7 


0.276 


0.011 



5.4.1. Erroneous distance 

A wrong assumption about the distance gives in a verti- 
cal shift in the CMD (all other parameters give a diago- 
nal shift) and cannot be masked out through a correlated 
change of any of the other parameters. This results in quite 
a wide range in the values of the other parameters, except 
the value of a. The unexpected result is that irrespective 
if the distance is too short or too far: the power-law IMF 
slope flattens! 

If the distance is overestimated there are more synthetic 
stars present at fainter magnitudes. To get relatively more 
synthetic stars at brighter magnitudes one needs to flatten 
the power-law IMF slope. 

If the distance is underestimated then more synthetic stars 
are present at brighter magnitudes. One expects that a 
steeper power-law IMF slope is required as compensation. 



lower end of the main sequence. They are taken away from 
the stars located at brighter and brighter magnitudes. One 
therefore requires also in this case a flatter IMF slope. 

A flatter slope of the power-law IMF can be a hint that 
the distance of the stellar aggregate is wrong. Or it might 
be a hint that the zero point of the adopted synthetic pho- 
tometric system is different from the actual photometric 
system used. 

Recognizing that the slope is indeed flatter than the 
majority of the other cases outlined in Table ^ one may 
start to explore the assumption that the distance is wrong: 
release the constraint during the next exploration. 

The sensitivity to the distance implies that AMORE 
can be used to determine the distance to a stellar aggre- 
gate quite reliably. A bonus is that due to an initially 
wrongly assumed distance the extinction is in most cases 
better constrained. 

5.4.2. Erroneous extinction, age and metallicity 

In (V,V-I) CMDs a strong correlation between e xtinct ion, 
age and metallicity exists (see also Ng & Bertelli 1996 and 
references cited therein). A higher value of the extinction 
can be compensated by a younger age and/or a lower met- 
allicity. Indeed, the results in Tabic A. 3 show that this ac- 
tually occurs for the lower age and metallicity limit. The 
upper age and metallicity limit, however, drifts away in 
the opposite direction to compensate for 'erroneous cor- 
rections' applied by other parameters. 

A higher value of the star formation index results in a 
lower number of stars at the upper age metallicity limit. 
To get a sufficient number of high metallicity stars one has 
to stretch the upper metallicity limit to a slightly higher 
value. 

We further notice that a wrong value for the age and 
metallicity does not affect the extinction significantly. Our 
findings indirectly supports the method to determine high 
resolution (4' x 4') extinction maps towards the Galactic 
bulge by Schultheis et al. ( |1999 ) with the data obtained 
for the DeNIS project (Epchtcin et al. |1997D . 



This is not always true, see Table A. 3. Stars pop up at the 
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5.4.3. Erroneous slope for power-law IMF 

The slope of the power-law IMF is very strongly con- 
strained (Ng 1998) for the test population. As a conse- 
quence the changes in the values of the remaining param- 
eters are not extremely large. The slightly larger value 
of the slope a pushes slightly more synthetic stars to 
fainter magnitudes, introducing a relative deficiency of 
stars at brighter magnitudes. This is compensated through 
a younger age, a decrease of the lower metallicity limit and 
an increase of the upper metallicity limit. The higher value 
for the upper age and metallicity limit compensates in its 
turn for the overestimation of the exponential star forma- 
tion index. 

A lower value of a is partly compensated for by lowering 
the upper age limit, lowering the lower metallicity limit 
and overestimating the upper metallicity limit. 



5.4.4. Erroneous index for exponential SFR 

A larger index for an exponentially decreasing SFR pushes 
more stars of the population to the blue edge of the CMD, 
resulting in a slightly bluer stellar population. AMORE 
compensates this by mainly increasing the upper metalli- 
city limit, i.e. reddening the synthetic stellar population. 
A lower value for the SFR index has the opposite effect. 
AMORE compensates for the now slightly redder popula- 
tion by lowering the upper metallicity limit, making the 
population bluer. 



6. Discussion 

6.1. Relative contribution of the parameters to the 
fitness 

Figure | and the Tables |, |, 0, |Aj & [Aj provide the 
hint that not all parameters have an equal contribution 
to the fitness. It appears that the f3 and Zhigh parame- 
ter can vary considerably and still yield a decent fitness. 
Moreover, Tables |A.2| & |5| indicate that knowledge of the 
value of the Zhigh parameter results in acceptable values 
for the other parameters. 

The origin of this behaviour lies in the implicit definition 
of the exponential star formation rate (for (3 = 1 one has 
a decreasing star formation towards a younger age) at- 
tached to a linear age-metallicity relation. The latter re- 
lation will give less metal-richer stars. The small number 
of stars with higher metallicity induces a larger variation 
of the Zhigh parameter without affecting significantly the 
overall fitness. 



6.2. Convergence 

The fine-tuning of the genetic algorithm is a tedious task. 
It is not straightforward to find the optimum setting for 
the problem to be solved. One has to balance the explor- 
ing quality through crossovers against the variation of the 
parameters through (creep) mutations. 



We did not want to deal with a mutation dominated 
search, because it tends to move farther away from an 
optimum parameter setting in the majority of the cases. 
We used therefore a relatively high crossover probability 
(pcross) and we set the mutations at a fixed rate, such 
that on average only 2.8 mutations occur in the gene pool 
of each individual. 

At a certain stage however one requires the varia- 
tion of other correlated parameters to obtain an impro- 
vement. This becomes particularly necessary when ap- 
proaching the optimum setting of the parameters. A 
favourable crossover and mutation might do the trick, 
but it might t ake a while before this occurs. We intro- 
duced in Sect. 2.6.5 the possibility that two parameters 
might be more sensitive to mutations than others. This 
approach gave better results for the majority of the trial 
cases (see Table ||) , but it failed to obtain improvements 
when changes of one parameter were neutralized through 
the variation of one or more parameters. The distance- 
extinction and the age-metallicity degeneracies slow down 
the convergence of AMORE for /> 0.3, see Fig. f|. 

One of the modifications to consider for future imple- 
mentation is a two-chromosome approach. In that case 
acceptable values for the parameters do not shift out of 
the population if the overall fitness is less, but still re- 
side in the gene pool as a recessive quality. This however, 
will require a major extension to PIKAIA and a significant 
amount of genetic research to be done about dominant 
and recessive qualities in the AMORE gene pool. 
Another modification to consider in order to improve the 
accuracy and to speed up convergence, is to replace the 
finite resolution of the digital encoding scheme with a ge- 
netic coding based on floating point, i.e. each gene on the 
chromosome is represented by one floating point number. 
According to Michalewicz ( 1996 ) a real encoding scheme 
can be superior and improve convergence. Such an encod- 
ing scheme is indeed to be included in the next release of 
PIKAIA 2.0 (Charbonneau; in preparation). 

6.2.1. Unstable solutions 

In one test (fixing Ay at one sigma above the original 
value for model 9) no convergence was achieved and the 
run was aborted. Because AMORE is quite sensitive to 
rounding these effects can be circumvented by slightly 
altering the input parameters. We decided in this case 
against such an action, because that would make the sam- 
ple inhomogeneous. 



6.3. Degeneracy 

Isochrones for a particular age and metallicity can be mim- 
icked with another set of isochrones of different age and 
metallicity (Worthey |1994| , Chariot et al. |1996[ , and ref- 
erences cited therein). This degeneracy of the parameter 
space increases if one considers the distance and the ex- 
tinction towards a stellar aggregate. There is no straight 
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Table 8. Description of the diagnostic statistics for simulations with setup parameters from model |A.l[ -40. One of the 
parameters is forced to a value ± la from its original value. See Sects. 5.4 and 2.1 for additional details 



model 


offset 


No, not 


Ns,not 


N match 


F x 


F P 


F 


fA 


ideal 












5000 


0.000 


0.000 


0.000 


1.000 


free 






60 


60 


4940 


0.847 


0.849 


1.438 


0.410 


fixed i 


log a i^pcj 


-Iff 


91 


91 


4909 


0.788 


1.287 


2.278 


0.305 






+ lff 


85 


85 


4915 


0.722 


1.202 


1.967 


0.337 


nxed: 


A , - 


+lcr 


98 


98 


4902 


0.932 


1.386 


3.790 


0.264 


fixed i 


log zi ow yyT) 


-Iff 


91 


91 


4901 


0.897 


1.287 


2.460 


0.289 






+ lff 


94 


94 


4906 


0.905 


1.329 


2.587 


0.278 


fixed: 


logth<,h(yr) 


-Iff 


106 


106 


4894 


0.977 


1.499 


3.202 


0.238 






+ 1(7 


126 


126 


4874 


1.017 


1.782 


4.211 


0.192 


fixed: 




-la 


78 


78 


4922 


0.907 


1.103 


2.043 


0.329 






+ 1(7 


80 


80 


4920 


0.835 


1.131 


1.977 


0.336 


fixed: 


[Z]hgh 


-la 


148 


148 


4852 


0.896 


2.094 


5.188 


0.162 






+la 


94 


94 


4906 


0.867 


1.329 


2.519 


0.284 


fixed: 


a 


-la 


95 


95 


4905 


0.886 


1.344 


2.590 


0.279 






+la 


103 


103 


4897 


0.894 


1.457 


2.921 


0.255 


fixed: 


P 


-la 


145 


145 


4855 


0.930 


2.051 


5.071 


0.165 






+la 


103 


103 


4897 


0.834 


1.457 


2.818 


0.262 



forward method to circumvent partial degeneracy of the 
parameters to be explored. One might consider to apply 
AMORE for the analysis of colour-colour diagrams in order 
to rule out the distance, to determine the extinction & a 
number of other parameters, and finally to determine the 
distance to the stellar population from one of the CMDs. 
The combined analysis of colour-magnitude and colour- 
colour diagrams is expected to improve the results ob- 
tained by AMORE so far. However, this requires that at 
least one additional colour should be available for each 
star considered above a certain detection and complete 
ness threshold. Moreover, as was mentioned in Sect. |5.3.1 



knowledge of the extinction does not automatically imply 
that the distance can be retrieved accurately. 

As demonstrated in Sect. [5.1.l| the degeneracy among 
parameters becomes noticeable for />0.25 or _F<3. This 
corresponds to systematic offset for each parameter of on 
average ~ 0.6(7^ and at maximum ~ lofc. The Poisson un- 
certainty of the original population results in a fitness of 
/ = 0.43 (F~1.33). However, solutions with a comparable 
fitness do exist due to the degeneracy of the parameter 
space. A direct consequence is that there is an intrinsic 
offset present among the parameters amounting to on av- 
erage ~ 0.4(7fc and at maximum ~ 0.7<7fc. This intrinsic 
offset is present in the solutions obtained with AMORE 
and actually is responsible for slowing down the conver- 
gence in the fitness range 0.30 </< 0.43. ft will therefore 
be nearly impossible to recover in one pass the original 
input values. However, some improvements might be ob- 
tained by averaging the parameter values obtained from 
AMORE runs with different initial conditions. 

7. Conclusions 

We demonstrate that an automatic search can be made for 
the astrophysical parameters of a synthetic stellar pop- 



ulation from the analysis of colour-magnitude diagrams 
with an optimizer, based on a genetic algorithm. However, 
AMORE tends to slightly underestimate the distance. It 
subsequently attempts to compensate this with an higher 
extinction, a higher upper age and a slightly steeper slope 
for the power-law IMF. At / > 0.3 the combined ef- 
fect of the age-mctallicity & the distance-extinction de- 
generacy slows down the convergence. The data suggests 
that AMORE has more problems dealing with the age- 
metallicity than with the distance-extinction degeneracy. 
For general purpose, however, the retrieved values are suf- 
ficiently accurate. 

8. Future work 

The good results obtained so far for a single synthetic 
stellar population is an indication about AMORE's poten- 
tial for the detailed analysis of CMDs. The next step is 
to improve one step at a time various aspects of AMORE 
before it can be used as an interpretative tool for large 
photometric surveys. 

Despite limitations in the input physics of the underlying 
stellar evolutionary tracks and by the transformation from 
the theoretical to the observational plane, the results with 
real data from Gallart et al. (1999), who uses the same set 



of evolutionary tracks, are encouraging. It will therefore 
be important to verify first with, for example, well studied 



open clusters (see Carraro et al. 1998, 199E and references 
cited therein) for which age and metallicity range we may 
apply AMORE safely. Extinction is also of some concern, 
because a high extinction may result in the MS turnoff 
point to fall below the detection limit. This would deprive 
AMORE of a clear reference point. 

Another case of interest is of course the question how 
many different stellar populations can be distinguished 
with AMORE. Separating multiple, mixed populations 



Y.K. Ng et al.: Automatic Observation Rendering (AMORE) 



19 



from each other through the automated and the objective 
analysis of colour-magnitude diagrams could be a valuable 
tool for the analysis of galaxy formation and evolution. 
This requires a rigorous follow up study on the separation 
of multiple (synthetic) populations. One further ought to 
verify if the automated analysis of colour-colour diagram 
can reduce the effects of error cancellation between dis- 
tance and extinction. 

Finally, after a succesful implementation, testing and 
validation phase, we plan to combine AMORE with the 



Padova spectrophotometric code (see Bressan et al. 1994 



1996; Tantalo et al. 1996, 1998ab). A synthetic population 



has to be generated, containing sub-populations with dif- 
ferent ages and mctallicitics. Then a synthetic spectrum 
must be generated for the mixed population and subse- 
quently used as input for a synthetic, spectral fitting pro- 
gram to determine the underlying stellar populations. In 
this way one can establish the calibration of the spec- 
trophotometric tool in a self-consistent way. Furthermore, 
an implicit verification can be made that the populations 
are consistent with those obtained from a CMD analysis 
with AMORE. 

Acknowledgements. The research was partly supported by 
TMR grant ERBFMRX-CT96-0086 from the European 
Community (Network: Formation and Evolution of Galaxies), 
by the Italian Ministry of University, Scientific Research and 
Technology (MURST) and by the Italian Space Agency (ASI). 

References 

Aparicio A., 1999, Proceedings IAU symposium 192 'The stel- 
lar content of the Local Group galaxies', 7-11 Sept. 1998, 
P. Whitelock and R. Cannon (eds.), ASP conference Series, 
304 

Banzhaf W., Nordin P., Keller R.E., Francone F.D., 1998, 

'Genetic Programming', (Morgan Kaufmann Publishers) 
Beichman C.A., Chester T.J., Skrutskie M., Low F.J., Gillett 

F, 1998, PASP 110, 480 
Bobinger A., 2000, A&A 357, 1170 
Bertelli G., Bressan A., Chiosi C, 1985, A&A 150, 33 
Bertelli G., Betto R., Bressan A., Chiosi C, Nasi E., Vallenari 

A., 1990, A&AS 85, 845 
Bertelli G., Mateo M., Chiosi C, Bressan A., 1992, ApJ 387, 

320 

Bertelli G., Bressan A., Chiosi C, Fagotto F., Nasi E., 1994, 

A&AS 106, 275 
Bertelli G., Bressan A., Chiosi C, Ng Y.K., Ortolani S., 1995, 

A&A 301, 381 

Bertelli G., Bressan A., Chiosi C, Ng Y.K., 1996, A&A 310, 
115 

Bessell M.S., Brett J.M., 1988, PASP 100, 1134 
Bressan A., Fagotto F., Bertelli G., Chiosi C, 1993, A&AS 
100, 647 

Bressan A., C. Chiosi, Fagotto F., 1994, ApJS 94, 63 
Bressan A., C. Chiosi, Tantalo R., 1996, A&A 311, 425 
Carraro G., Ng Y.K., Portinari L., 1998, MNRAS 296, 1045 
Carraro G., Girardi L., Chiosi C, 1999, MNRAS 309, 430 
Charbonneau P., 1995, ApJS 101, 309 

Charbonneau P., Knapp B., 1996, "A User's Guide to 
PIKAIA 1.0", NCAR Technical Note 418+IA (Boulder: 
National Center for Atmospheric Research) 



Charbonneau P., Tomczyk S., Shou J., Thompson M.J., 1998, 

ApJ 496, 1015 
Chariot S., Worthey G., Bressan A., 1996, ApJ 457, 625 
Chiosi C, 1996, in Proceedings 'From Stars to Galaxies', 9-13 

October, 1995, Elounda (Crete; Greece), C. Leitherer, U.F. 

von Alvensleben and J. Huchra (eds.), ASP Conference 

Series Vol. 98, 181 
Chiosi C, Bertelli G., Meylan G., Ortolani S., 1989, A&A 219, 

167 

da Costa L., 1997, Messenger 88, 34 

da Costa L., et al., 1998, Messenger 91, 49 

Davis L., 1991, Handbook of Genetic Algorithms (New York: 

Van Nostrand Reinhold) 
De Jong K.A., 1993, in 'Foundations of Genetic Algorithms 2', 

L.D. Whitley (ed.) (San Mateo; Morgan Kaufmann), 5 
Dolphin A., 1997, New Astronomy 2, 397 

Dolphin A.E., Walker A.R., Hodge P.W., et al, 2001, ApJ 562, 
303 

Dolphin A., 2002, MNRAS 332, 91 

Epchtein N., Batz de B., Capoani L., et al., 1997, ESO 

Messenger 87, 27 
Epchtein N., Deul E., Derriere S., et al., 1999, A&A 349, 236 
Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994a, A&AS 

104, 365 

Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994b, A&AS 

105, 29 

Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994c, A&AS 

105, 39 
Fan X., 1999, AJ 117, 2528 

Gallart C, Aparicio A., Bertelli G., Chiosi C, 1996, AJ 112, 
1950 

Gallart C, 1998, ApJ 495, L43 

Gallart C, Freedman W.L., Aparicio A., Bertelli G., Chiosi C, 

1999, AJ 118, 2245 
Geha M.C., Holtzman J.S., Mould J.R., et al., 1998, AJ 115, 

1045 

Gibson S.E., Charbonneau P., 1998, JGR 103, No. A7, 14511 
Gilmore G.F., Perryman M.A., Lindegren L., 1998, SPIE, 3350, 
541 

Goldberg D.E., 1989, Genetic Algorithms in Search, 
Optimization & Machine Learning, (Reading, MA: 
Addison- Wesley) 

Harris J., Zaritsky D., 2001, ApJS 136, 25 

Hernandez X., Valls-Gabaud D., Gilmore G., 1999, MNRAS 
304, 705 

Holland J.H., 1975, Adaptation in Natural and Artificial 
Systems (I s ed.; Ann Harbor: University of Michigan Press 
- 2 nd ed.; 1992, Cambridge: MIT Press) 

Holtzman J.A., Mould J.R., Gallagher J.S., et al., 1997, AJ 
113, 656 

Holtzman J.A., Gallagher J.S., Cole A.A., et al., 1999, AJ 118, 
2262 

Iglesias C.A., Rogers F.J., Wilson B.C., 1992, ApJ 397, 717 
Kaastra J.S., Mewe R., Liedahl D.A., et al., 1996, A&A 314, 
547 

Kennelly E.J., Walker G.A.H., Catala C, et al., 1996, A&A 
313, 571 

Lamontagne R., Moffat A.F.J., Drissen L., Robert C, 

Matthews J.M., 1996, AJ 112, 2227 
Mcintosh S.W., Diver D.A., Judge P.G., et al., 1998, A&AS 

132, 145 

Metcalfe T.S., 1999, AJ 117, 2503 

Mewe R., Kaastra J.S., White S.M., Pallavich R., 1996, A&A 
315, 170 



20 



Y.K. Ng et al.: Automatic Observation Rendering (AMORE) 



Michalewicz Z., 1996, 'Genetic algorithms + data structures fourth test, in which one parameter was fixed one sigma 

= evolution programs', 3 rd revised and extended edition from its correct value as discussed in Sect. 4.3.4, arc given 

(Springer- Verlag) m Table [Aj| 

Ng Y.K., 1994, Ph.D. thesis, Leiden University, the The tables in the Appendix are available in electronic form 



Netherlands 

Ng Y.K., 1997a, in Proceedings 'The impact of large-scale near- 
IR sky surveys', 24-26 April, 1996, Puerto de la Cruz 
(Tenerife; Spain), P. Garzon-Lopez et al. (eds.), 91 

Ng Y.K., 1997b, in Proceedings 12th IAP Colloquium 'Variable 
stars and the astrophysical returns from microlensing sur- 
veys', 8-12 July, 1996, Paris (France), R. Ferlet, J.-P. 
Maillard and B. Raban (eds.), 113 

Ng Y.K., 1998, A&AS 132, 133 

Ng Y.K., Bertelli G., 1996, A&A 315, 116 

Ng Y.K., Bertelli G., Bressan A., Chiosi C, Lub J., 1995, A&A 
295, 655 (erratum A&A 301, 318) 

Ng Y.K., Bertelli G., Chiosi C, Bressan A., 1996, A&A 310, 
771 

Ng Y.K., Bertelli G., Chiosi C, Bressan A., 1997, A&A 324, 
65 

Noyes R.W., Lha S., Korzennik S.G., et al, 1997, ApJ 483, 
Llll 

Olsen K.A.G., 1999, AJ 117, 2244 

Paczyhski B., Udalski A., Szymanski M., 1999, Acta Astron. 
49, 319 

Perryman M.A.C., de Boer K.S., Gilmore G., et al., 2001, A&A 
369, 339 

Pagel B.E.J. , Portinari L., 1998, MNRAS 298, 747 
Powell M.J.D., 1964, Comput.J. 7, 155 

Press W.H., Flannery B.P., Teukolsky S.A., Vettering 
W.T., 1986, Numerical Recipes (2nd edition), Cambridge 
University Press 

Renzini A., da Costa L., 1997, Messenger 87, 23 

Rieke G.H., Lebofsky M.J., 1985, ApJ 288, 618 

Saha P., 1998, ApJ 115, 1206 

Schultheis M., Ganesh S., Simon G., et al, 1999, A&A 349, 
L69 

Skrutskie M., 1998, in 'The impact of near-infrared sur- 
veys on galactic and extragalactic astronomy', Proc. 3 rd 
Euroconference, Kluwer ASSL 230, 11 

Stetson P.B., Harris W.E., 1988, AJ 96, 909 

Tantalo R., Chiosi C, Bressan A., Fagotto F., 1996, A&A 311, 
361 

Tantalo R., Chiosi C, Bressan A., 1998, A&A 333, 419 
Tantalo R., Chiosi C, Bressan A., Marigo P., Portinari L., 

1998, A&A 335, 823 
Udalski A., Kubiak M., Szymanski M., 1997, Acta Astron. 47, 

319 

van de Hulst H.C., 1949, Rech.Astr.Obs.Utrecht, Vol. 11 part 
2 

Wahde M., 1998, A&AS 132, 417 
Worthey G., 1994, ApJS 95, 107 



at the CDS. 



Appendix A: AMORE test results, the data 



The appendix contains the data of the simulations as dis- 
cussed in Sect. f|. The first test, in which we explored 
different values for the PIKAIA control parameters, is de- 



scribed in Sect. 4.3.1 and the data are given in Table A.l 



The data of the third test, in which one of the astrophys- 
ical parameters was fixed at its correct value as discussed 
in Sect. 4.3.3, are given in Table A. 2. The data on the 
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2.6.3| , |2.6.4 |2.6.5| and |J for 



Tabic A.l: AMORE with the test population as specified in Sect. 4.3.1 . See 2.6.2 , 
additional details, pcross is the crossover probability, rcross is the multi-point crossover rate, rbrood determines 
the amount of offspring two parents can have, pcreep is the creep mutation rate and pcorr the correlated mutation 
rate. The generation number indicates at which generation the resulting fitness value /a obtained with AMORE first 
emerged 



model 


pcross 


rcross 


rbrood 


pcreep 


pcorr 


fA 


generation 


1 


0.85 


1.00 


1 


00 


0.0 


0.0 


0.17266 


361 


2 


0.85 


1.00 


1 


00 


0.3 


0.0 


0.17266 


361 


3 


0.85 


1.00 


1 


00 


0.7 


0.0 


0.17266 


361 


4 


0.85 


1.00 


1 


00 


0.0 


0.3 


0.27551 


221 


5 


0.85 


1.00 


1 


00 


0.0 


0.7 


0.26830 


221 


6 


0.85 


1.00 


1 


00 


0.3 


0.3 


0.30488 


381 


7 


0.85 


1.00 


1 


00 


0.7 


0.3 


0.32005 


341 


8 


0.85 


1.00 


1 


00 


0.3 


0.7 


0.32167 


201 


9 


0.85 


1.00 


1 


00 


0.7 


0.7 


0.38310 


341 


10 


0.50 


1.00 


1 


00 


0.0 


0.0 


0.35662 


381 


11 


0.50 


1.00 


1 


00 


0.3 


0.0 


0.35662 


381 


12 


0.50 


1.00 


1 


00 


0.7 


0.0 


0.35662 


381 


13 


0.50 


1.00 


1 


00 


0.0 


0.3 


0.34961 


281 


14 


0.50 


1.00 


1 


00 


0.0 


0.7 


0.17008 


261 


15 


0.50 


1.00 


1 


00 


0.3 


0.3 


0.36025 


381 


16 


0.50 


1.00 


1 


00 


0.7 


0.3 


0.24048 


261 


17 


0.50 


1.00 


1 


00 


0.3 


0.7 


0.23956 


321 


18 


0.50 


1.00 


1 


00 


0.7 


0.7 


0.29463 


321 


19 


0.85 


2.00 


1 


00 


0.0 


0.0 


0.31700 


384 


20 


0.85 


2.00 


1 


00 


0.3 


0.0 


0.31700 


384 


21 


0.85 


2.00 


1 


00 


0.7 


0.0 


0.31700 


384 


22 


0.85 


2.00 


1 


00 


0.0 


0.3 


0.25740 


341 


23 


0.85 


2.00 


1 


00 


0.0 


0.7 


0.35805 


341 


24 


0.85 


2.00 


1 


00 


0.3 


0.3 


0.34120 


181 


25 


0.85 


2.00 


1 


00 


0.7 


0.3 


0.31290 


141 


26 


0.85 


2.00 


1 


00 


0.3 


0.7 


0.30913 


341 


27 


0.85 


2.00 


1 


00 


0.7 


0.7 


0.32212 


261 


28 


0.50 


2.00 


1 


00 


0.0 


0.0 


0.16122 


261 


29 


0.50 


2.00 


1 


00 


0.3 


0.0 


0.16122 


261 


30 


0.50 


2.00 


1 


00 


0.7 


0.0 


0.16122 


261 


31 


0.50 


2.00 


1 


00 


0.0 


0.3 


0.36796 


381 


32 


0.50 


2.00 


1 


00 


0.0 


0.7 


0.28494 


381 


33 


0.50 


2.00 


1 


00 


0.3 


0.3 


0.27263 


381 


34 


0.50 


2.00 


1 


00 


0.7 


0.3 


0.16352 


301 


35 


0.50 


2.00 


1 


00 


0.3 


0.7 


0.35149 


301 


36 


0.50 


2.00 


1 


00 


0.7 


0.7 


0.34031 


381 


37 


0.85 


3.00 


1 


00 


0.0 


0.0 


0.38297 


301 


38 


0.85 


3.00 


1 


00 


0.3 


0.0 


0.38297 


301 


39 


0.85 


3.00 


1 


00 


0.7 


0.0 


0.38297 


301 


40 


0.85 


3.00 


1 


00 


0.0 


0.3 


0.41019 


341 


41 


0.85 


3.00 


1 


00 


0.0 


0.7 


0.32012 


381 


42 


0.85 


3.00 


1 


00 


0.3 


0.3 


0.36840 


341 


43 


0.85 


3.00 


1 


00 


0.7 


0.3 


0.31281 


361 


44 


0.85 


3.00 


1 


00 


0.7 


0.3 


0.31281 


361 


45 


0.85 


3.00 


1 


00 


0.3 


0.7 


0.30571 


381 


46 


0.50 


3.00 


1 


00 


0.0 


0.0 


0.15614 


301 


47 


0.50 


3.00 


1 


00 


0.3 


0.0 


0.15614 


301 


48 


0.50 


3.00 


1 


00 


0.7 


0.0 


0.15614 


301 


49 


0.50 


3.00 


1 


00 


0.0 


0.3 


0.23392 


301 


50 


0.50 


3.00 


1 


00 


0.0 


0.7 


0.16274 


221 


51 


0.50 


3.00 


1 


00 


0.3 


0.3 


0.25123 


281 
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model 


pcross 


rcross 


rbrood 


pcreep 


pcorr 


fA 


generation 


52 


0.50 


3.00 


1.00 


0.7 


0.3 


0.27151 


261 


53 


0.50 


3.00 


1.00 


0.3 


0.7 


0.35597 


301 


54 


0.50 


3.00 


1.00 


0.7 


0.7 


0.27743 


361 


55 


0.85 


1.00 


2.00 


0.0 


0.0 


0.36221 


161 


56 


0.85 


1.00 


2.00 


0.3 


0.0 


0.36221 


161 


57 


0.85 


1.00 


2.00 


0.7 


0.0 


0.36221 


161 


58 


0.85 


1.00 


2.00 


0.0 


0.3 


0.16284 


261 


59 


0.85 


1.00 


2.00 


0.0 


0.7 


0.30251 


381 


60 


0.85 


1.00 


2.00 


0.3 


0.3 


0.32133 


361 


61 


0.85 


1.00 


2.00 


0.7 


0.3 


0.33565 


361 


62 


0.85 


1.00 


2.00 


0.3 


0.7 


0.27452 


341 


63 


0.85 


1.00 


2.00 


0.7 


0.7 


0.34660 


301 


64 


0.50 


1.00 


2.00 


0.0 


0.0 


0.27683 


381 


65 


0.50 


1.00 


2.00 


0.3 


0.0 


0.27683 


381 


66 


0.50 


1.00 


2.00 


0.7 


0.0 


0.27683 


381 


67 


0.50 


1.00 


2.00 


0.0 


0.3 


0.16833 


201 


68 


0.50 


1.00 


2.00 


0.0 


0.7 


0.38586 


321 


69 


0.50 


1.00 


2.00 


0.3 


0.3 


0.34390 


361 


70 


0.50 


1.00 


2.00 


0.7 


0.3 


0.33531 


241 


71 


0.50 


1.00 


2.00 


0.3 


0.7 


0.32646 


221 


72 


0.50 


1.00 


2.00 


0.7 


0.7 


0.34591 


381 


73 


0.85 


1.00 


4.00 


0.0 


0.0 


0.34607 


221 


74 


0.85 


1.00 


4.00 


0.3 


0.0 


0.34607 


221 


75 


0.85 


1.00 


4.00 


0.7 


0.0 


0.34607 


221 


76 


0.85 


1.00 


4.00 


0.0 


0.3 


0.32143 


361 


77 


0.85 


1.00 


4.00 


0.0 


0.7 


0.27159 


381 


78 


0.85 


1.00 


4.00 


0.3 


0.3 


0.26001 


221 


79 


0.85 


1.00 


4.00 


0.7 


0.3 


0.19802 


241 


80 


0.85 


1.00 


4.00 


0.3 


0.7 


0.33032 


381 


81 


0.85 


1.00 


4.00 


0.7 


0.7 


0.17799 


341 


82 


0.50 


1.00 


4.00 


0.0 


0.0 


0.27728 


201 


83 


0.50 


1.00 


4.00 


0.3 


0.0 


0.27728 


201 


84 


0.50 


1.00 


4.00 


0.7 


0.0 


0.27728 


201 


85 


0.50 


1.00 


4.00 


0.0 


0.3 


0.35783 


381 


86 


0.50 


1.00 


4.00 


0.0 


0.7 


0.28070 


341 


87 


0.50 


1.00 


4.00 


0.3 


0.3 


0.33717 


381 


88 


0.50 


1.00 


4.00 


0.7 


0.3 


0.18275 


341 


89 


0.50 


1.00 


4.00 


0.3 


0.7 


0.24795 


361 


90 


0.50 


1.00 


4.00 


0.7 


0.7 


0.28384 


341 


91 


0.85 


2.00 


2.00 


0.0 


0.0 


0.29274 


261 


92 


0.85 


2.00 


2.00 


0.3 


0.0 


0.29274 


261 


93 


0.85 


2.00 


2.00 


0.7 


0.0 


0.29274 


261 


94 


0.85 


2.00 


2.00 


0.0 


0.3 


0.24259 


381 


95 


0.85 


2.00 


2.00 


0.0 


0.7 


0.26513 


161 


96 


0.85 


2.00 


2.00 


0.3 


0.3 


0.33411 


361 


97 


0.85 


2.00 


2.00 


0.7 


0.3 


0.25808 


121 


98 


0.85 


2.00 


2.00 


0.3 


0.7 


0.35686 


341 


99 


0.85 


2.00 


2.00 


0.7 


0.7 


0.37254 


399 


100 


0.50 


2.00 


2.00 


0.0 


0.0 


0.28259 


361 


101 


0.50 


2.00 


2.00 


0.3 


0.0 


0.28259 


361 


102 


0.50 


2.00 


2.00 


0.7 


0.0 


0.28259 


361 


103 


0.50 


2.00 


2.00 


0.0 


0.3 


0.34559 


341 


104 


0.50 


2.00 


2.00 


0.0 


0.7 


0.39464 


320 


105 


0.50 


2.00 


2.00 


0.3 


0.3 


0.27112 


181 


106 


0.50 


2.00 


2.00 


0.7 


0.3 


0.29866 


321 
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model 


pcross 


rcross 


rbrood 


pcreep 


pcorr 


fA 


generation 


107 


0.50 


2.00 


2.00 


0.3 


0.7 


0.23660 


181 


108 


0.50 


2.00 


2.00 


0.7 


0.7 


0.35769 


261 


109 


0.85 


2.00 


4.00 


0.0 


0.0 


0.32628 


181 


110 


0.85 


2.00 


4.00 


0.3 


0.0 


0.32628 


181 


111 


0.85 


2.00 


4.00 


0.7 


0.0 


0.32628 


181 


112 


0.85 


2.00 


4.00 


0.0 


0.3 


0.26910 


361 


113 


0.85 


2.00 


4.00 


0.0 


0.7 


0.37278 


389 


114 


0.85 


2.00 


4.00 


0.3 


0.3 


0.31434 


241 


115 


0.85 


2.00 


4.00 


0.7 


0.3 


0.18222 


321 


116 


0.85 


2.00 


4.00 


0.3 


0.7 


0.26176 


341 


117 


0.85 


2.00 


4.00 


0.7 


0.7 


0.19495 


241 


118 


0.50 


2.00 


4.00 


0.0 


0.0 


0.17656 


361 


119 


0.50 


2.00 


4.00 


0.3 


0.0 


0.17656 


361 


120 


0.50 


2.00 


4.00 


0.7 


0.0 


0.17656 


361 


121 


0.50 


2.00 


4.00 


0.0 


0.3 


0.17248 


181 


122 


0.50 


2.00 


4.00 


0.0 


0.7 


0.27790 


381 


123 


0.50 


2.00 


4.00 


0.3 


0.3 


0.17117 


221 


124 


0.50 


2.00 


4.00 


0.7 


0.3 


0.36171 


381 


125 


0.50 


2.00 


4.00 


0.3 


0.7 


0.38431 


392 


126 


0.50 


2.00 


4.00 


0.7 


0.7 


0.31711 


400 


127 


0.85 


3.00 


2.00 


0.0 


0.0 


0.16288 


161 


128 


0.85 


3.00 


2.00 


0.3 


0.0 


0.16288 


161 


129 


0.85 


3.00 


2.00 


0.7 


0.0 


0.16288 


161 


130 


0.85 


3.00 


2.00 


0.0 


0.3 


0.34149 


341 


131 


0.85 


3.00 


2.00 


0.0 


0.7 


0.33630 


281 


132 


0.85 


3.00 


2.00 


0.3 


0.3 


0.24355 


301 


133 


0.85 


3.00 


2.00 


0.7 


0.3 


0.33731 


221 


134 
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rcross 


rbrood 
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fA 
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Table A. 2. Results of running AMORE with one of the astrophysical parameters fixed at its original value. These 
tests provide an indication about the influence of one parameter on the retrieval of the remaining parameters, see 
Sect. 4.3.3| and 5.3 for additional details 
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Table A. 3: Results of AMORE with one of the parameters fixed at one sigma from its original value. These tests provide 
an indication about the influence of one parameter on the retrieval of the remaining parameters. In p articu lar, how 
the remaining parameters balance this mis-match by moving away from their optimal value. See Sect, i .3.4 and 5.4 
for additional details 
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52 


0.32236 
0.19541 
0.19646 
n i 09m 

0.32867 
0.28939 


3.89424 
3.89205 
3.89803 

O.OiJ i oU 

3.89484 
3.89843 


0.060 
0.085 
0.068 
u.uoo 
0.038 
0.035 


9.92490 
9.80620 
9.80556 

O 7fi^S9 
y . / DOoi 

9.90919 
9.88752 


9.94949 
9.95304 
9.95382 

Q QA 81 fl 

9.95602 
9.95511 


-0.65235 
-0.65235 
-0.65235 

— U.OOZoO 

-0.65235 
-0.65235 


0.09935 
0.44787 
0.48659 

n A 8 7/1 9 

U.4:0 / 4:Z 

0.18180 
0.22274 


2.364 
2.367 
2.368 

O QIC 
Z.OlO 

2.356 
2.337 


0.023 
4.219 
4.525 
a no*? 

4.UZO 

0.601 
0.836 




+ la 


9 
14 
22 

i.A 

40 
52 


0.30963 
0.32883 
0.32293 

U.ZDO i O 

0.33593 
0.32988 


3.89302 
3.89295 
3.89110 
o.oyozo 
3.89840 
3.89341 


0.019 
0.022 
0.024 

U.UZ4 

0.032 
0.033 


9.89547 
9.91000 
9.91831 

O 8fi/l 1 A 

y .004:04: 

9.90057 
9.90827 


9.97009 
9.97315 
9.96435 
y.yoooy 
9.96348 
9.97162 


-0.55177 
-0.55177 
-0.55177 
n ^ 77 

-0.55177 
-0.55177 


0.26207 
0.19130 
0.18106 

n ^79^ ^ 
U . O (ZOO 

0.17853 
0.20534 


2.372 
2.365 
2.369 

Z.oOl 

2.349 
2.376 


1.833 
0.778 
0.671 
2 928 
0.885 
1.416 


log Zhigh 


-la 


9 
14 
22 

04 

40 
52 


0.29612 
0.29094 
0.27692 

V.ZoiioV 

0.16160 
0.30348 


3.89604 
3.89708 
3.89364 

Q OQ7QC 

3.90196 
3.89643 


0.077 
0.082 
0.066 
0.083 
0.043 
0.078 


9.92807 
9.92897 
9.94038 

o noenc 

y.yzBUo 
9.82928 
9.92740 


9.93582 

9.93400 

9.95069 

y.youo i 

10.10576 

9.93508 


-0.66338 
-0.67421 
-0.66032 
— u.Ooozy 
-0.79270 
-0.66677 


0.04239 
0.04239 
0.04239 

n n 1 9*? o 
U.U4zoy 

0.04238 
0.04239 


2.356 
2.359 
2.387 
2.357 
2.333 
2.355 


-0.232 
-0.139 
-0.337 
— 0.095 
-1.676 
-0.238 




+ la 


9 
14 
22 
34 
40 
52 


0.26702 
0.29357 
0.30635 
95823 
0.28417 
0.28883 


3.88982 
3.89392 
3.89366 
3 89389 
3.89886 
3.89501 


0.031 
0.012 
0.032 
0.046 
0.028 
0.020 


9.88500 
9.88192 
9.87176 
9 87585 

O . O J OOO 

9.85863 
9.88098 


9.97615 
9.97494 
9.96383 
9 96245 
9.95559 
9.96677 


-0.56217 
-0.53776 
-0.53473 
— 53400 
-0.57966 
-0.55136 


0.30979 
0.30979 
0.30979 
30979 
0.30979 
0.30979 


2.384 
2.371 
2.351 
2.357 
2.333 
2.367 


2.280 
2.129 
2.194 
2.578 
1.926 
2.284 


a 


-la 


9 
14 
22 
34 
40 
52 


0.24524 
0.28147 
0.27797 
0.22427 
0.27857 
0.28640 


3.90004 
3.90272 
3.89896 
3.89644 
3.90056 
3.90139 


0.053 
0.020 
0.042 
0.032 
0.042 
0.039 


9.85726 
9.87188 
9.86368 
9.82350 
9.86375 
9.85871 


9.94126 
9.94707 
9.94580 
9.94534 
9.94608 
9.94768 


-0.56638 
-0.54698 
-0.54982 
-0.50532 
-0.59271 
-0.58962 


0.28025 
0.24944 
0.26102 
0.45532 
0.26364 
0.27721 


2.316 
2.316 
2.316 
2.316 
2.316 
2.316 


1.932 
1.230 
1.511 
3.679 
1.294 
1.438 




+ la 


9 
14 
22 
34 
40 
52 


0.23946 
0.25586 
0.28357 
0.29934 
0.25505 
0.30437 


3.89295 
3.89172 
3.89109 
3.89146 
3.89153 
3.89128 


0.062 
0.014 
0.021 
0.023 
0.009 
0.026 


9.87606 
9.86873 
9.91357 
9.91497 
9.80529 
9.90426 


9.96018 
9.97973 
9.97110 
9.97387 
9.97968 
9.97836 


-0.57093 
-0.53518 
-0.56888 
-0.58765 
-0.52454 
-0.55141 


0.35127 
0.35320 
0.24524 
0.23086 
0.51555 
0.23875 


2.384 
2.384 
2.384 
2.384 
2.384 
2.384 


3.499 
2.919 
1.488 
1.461 
4.820 
1.620 



Y.K. Ng et al.: Automatic Observation Rendering (AMORE) 



Tabic A. 3: continued 



parameter 


offset 


model 


fA 


log d 


Av 


logtiom 


logth g /i 


[Z]lou, 




a 


P 


p 


-la 


9 


0.29110 


3.89717 


0.083 


9.93098 


9.93961 


-0.68243 


0.02188 


2.365 


-0.397 






14 


0.17266 


3.89314 


0.001 


9.83767 


10.08181 


-0.41912 


0.06629 


2.372 


-0.397 






22 


0.30211 


3.89619 


0.084 


9.91729 


9.94002 


-0.68522 


0.06730 


2.348 


-0.397 






34 


0.16366 


3.89491 


0.015 


9.81889 


10.08661 


-0.50001 


0.09684 


2.363 


-0.397 






40 


0.16472 


3.89974 


0.018 


9.80958 


10.07264 


-0.54158 


0.11992 


2.357 


-0.397 






52 


0.29684 


3.89592 


0.074 


9.91808 


9.93138 


-0.66681 


0.08493 


2.339 


-0.397 




+l<y 


9 


0.26953 


3.89638 


0.045 


9.86719 


9.95699 


-0.57054 


0.30738 


2.346 


2.397 






14 


0.29111 


3.89398 


0.023 


9.87948 


9.97371 


-0.54141 


0.30286 


2.371 


2.397 






22 


0.28842 


3.89548 


0.029 


9.86193 


9.95692 


-0.53941 


0.33240 


2.341 


2.397 






34 


0.27336 


3.89607 


0.048 


9.88719 


9.95851 


-0.56455 


0.27131 


2.355 


2.397 






40 


0.26195 


3.89766 


0.038 


9.87053 


9.94793 


-0.58042 


0.30908 


2.345 


2.397 






52 


0.27328 


3.89034 


0.022 


9.87427 


9.97430 


-0.54963 


0.33427 


2.375 


2.397 



